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The cross section of elastic electron-proton scattering taking place in an electron gas is calculated 
within the Closed Time Path method by means of the resummation of the one-loop self-energy in the 
photon propagator. Back-reaction of the colliding particles on the gas has been taken into account 
and is found to be dominating the cross section when the energy exchange falls in the vicinity of 
the Fermi energy. Back-reaction reflects the colliding particles-gas entanglement and makes the 
- - . colliding particle state mixed. The softness of the asymptotic particle-hole states of the gas makes 

■ the colliding particle trajectories more consistent and the collision irreversible, rendering in this 
I manner the scattering more classical in this regime. 
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I. INTRODUCTION 

Scattering experiments are imagined as idealized processes, described by transition probabilities between asymptotic 
in- and out-states without interactions, interpolated in time by turning on and off adiabatically interactions. In reality 
these processes take place in an environment that interacts with the system even if the final state interactions are 
JL . neglected. The influence of the environment on the asymptotic states can be taken into account in a simple manner 
as long as the environment possesses no asymptotic states. For instance, virtual vacuum polarizations might be 
, considered as environment from the point of view of the colliding particles and their influence on the final state can 
'— '■ be accounted by using dressed, quasifree particles in the asymptotic states. The influence of the environment on the 
cllision process is then restricted to the collision zone and can easily be recovered by the usual perturbation expansion. 
This problem becomes more interesting when the environment has asymptotic states. The difference between 

■ system and environment is that the former evolves from a fixed initial state to a given final one, the latter has initial 
^sj , conditions only: it is not identified by the detectors hence it follows an open-ended, undisturbed time evolution. 
T-H ■ One can, in principle, take such a possibility into account within the usual description of scattering, by summing 
C^l , the transition probability over possible orthogonal environment final states. But, this is not feasible in practice 

■ where this summation is carried out explicitly due to the large number of possible final environment states. A way 
T-H [ out from this problem is to build in the trace over environment states into the transition probability form the very 
^— I • beginning. The closed time path (CTP) method proposed by J. Schwinger [1] proves to be the appropriate framework 

I for this problem. It is not only tailored to calculate expectation value but transition probability can be recast in this 
^ . form. In addition, it contains the trace we need, over the Fock space of the theory to arrive at probabilities instead 
of transition amplitudes. An additional bonus of using this scheme is a simple way to build up the entanglement 
between the colliding particles and their environment and the mixing in the reduced density matrix for the scattering 

■ particles. 

■ " ' The goal of this paper is the description of the electron-proton elastic scattering taking place in a homogeneous 
electron gas. We consider the gas at vanishing temperature for the sake of simplicity; therefore, the model starts with 
a pure state, a Fermi-sphere of electrons, no photons, and an incident proton, at the initial time. The final state is 
identified by the momentum of the proton and a recoiled electron, being originally at rest with respect to the gas. 
This problem may be considered as the first step toward a more realistic description of collision processes where both 
the beam and the target are systems consisting of a large number of particles. Our main result is the form of the cross 
section as a sum of the expected expression containing the photon propagator dressed by the environment and another 
term, comprising the effects of colliding particles-environment entanglement. There are kinematical regions where this 
latter is more important than the former. We monitor the consistency [2-4] of the colliding particle trajectories and 
the dynamical breakdown of time reversal invariance, two necessary conditions for classical behavior. We find that 
both of them are strong for scattering where the entanglement contribution is dominant, supporting the expectation 
about the importance of entanglement in classical limit. 

The CTP method has already been used in different forms with different goals. It emerges in a thorough analysis of 
many-body theory [5]; it was introduced in solid state physics [6] to calculate retarded Green functions. The discussion 
of quantum field theory at finite temperature and/or finite density is based on this formalism [7-9] and algebraically 
defined thermal field theory (TFT) [10] can be derived from it, as well. Its path integral setting [11] opened the way 
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for a large number of applications [12]. Relaxation has been studied in an extensive manner [13] by deriving quantum 
kinetic equations in two different ways. One is to use the Schwingcr-Dyson, alias Kadanoff-Baym equations [5] to 
arrive at a Boltzmann-equation-like balance condition in a plasma [14-20]. Another way is to construct the analogy 
of a quantum Brownian motion [21], a Langevin-like equation in field theory [22-24]. 

The problem addressed here: to follow a single collision in an environment is a nonequilibrium process and requires 
the use of the CTP formalism that is equivalent with a special, unusual realization of TFT where one has access to 
transition probabilities and can trace the effects of entanglement between the colliding particles and their environment 
and estimate the consistency of the collision. The cross section is obtained by relying on the usual a simplification, 
adiabatic switching. It is assumed that interactions arc weak when the proton enters into the electron gas and when 
it leaves together with the kicked electron. Apart from some initial steps in TFT [25-27], we are not aware of the 
systematic application of the CTP scheme to describe transition probabilities or scattering cross sections. Another 
novel aspect of the results presented below is that the back- reaction of the collision process on the environment, a 
rather complicated effect that is usually ignored in the quantum kinetic approach, is retained. It can be taken into 
account in a simple manner by means of an explicit, algebraic solution of the Kadanoff-Baym equation. The resulting 
full-fledged interaction with balanced action-reaction forces leads to colliding charge-gas entanglement. Such an 
entanglement is incorporated in the photon self-energy, representing the coupling of photon to particle-hole excitations. 
The entanglement contributions dominate the cross section in the kinematical regime where the asymptotic excitations 
of the electron gas have sufficiently soft excitation spectrum. This softness makes entanglement nonperturbative and 
suppresses the quantum interference of trajectories thereby driving the collision classical. 

We start in Section II with the outline of an idealized version of an electron-proton scattering in an electron gas. 
The way of recognizing entanglement, mixing in the reduced density matrix and consistency is explained in Section 
III. The reduction formula and its perturbative evaluation is discussed in Section IV. The partially resummed photon 
propagator, appearing in the expression for the transition probability and its spectral functions are obtained in Section 
V. The numerical results for the cross section are presented in Section VI. Finally, the summary of our findings is 
given in Section VII. 

II. COLLISION IN A GAS 

We consider in this work elastic electron-proton scattering taking place in a homogeneous electron gas at vanishing 
temperature and finite density. A proton falls on the gas and electrons are knocked out from the gas by the collision. 
Though the intial state of electrons participating in the collision is not controllable we assume that those scattering 
events are kept where the initial electron was at rest with respect to the gas. 

The adiabatic switching hypothesis is more involved in this case than for scatterings in the vacuum because of the 
interacting environment. Let us distinguish the following well-separated times: ti <C tin ^ tout tfmai- First a 
noninteracting electron and proton beam are created at the initial time ti. The protons enter into the electron gas at 
tin , the colliding particles leave the gas at tout , and finally, the assembly of noninteracting final states of the colliding 
particles and the gas is removed at the final time t final- The adiabatic switching of the electromagnetic interaction 
should satisfy two constraints. On the one hand, the interactions should be strong enough to build up or deconstruct 
the interacting electron gas before tin or after tout, respectively. On the other hand, the interactions should be weak 
enough at tin and tout to render the entering into and leaving from the electron gas adiabatic, leaving behind no 
excitations. 

The initial state = ^p^aej^gas) written as a product of the operators 6],. and al., which create the initial 
proton and electron states, respectively and 

l*sa.)= n (1) 

denotes the state of the free gas with a hole at the quantum number e^. The amplitude of transition to a state 
containing a free proton and electron with quantum numbers pf and e/ is 

An = {n\aefbp^U{tf,ti)\^i). (2) 

where n denotes the quantum number of the remaining electron-photon component of the final state and U{tf,ti) 
stands for the time evolution operator. It is assumed that the detectors register the state of the colliding particles 
only and the electrons in the rest of the gas and photons are left in an arbitrary final state. The probability of the 
observed transition is therefore the sum over possible states of the remainder of the electron gas and photons, 

P{f^i) = J2\-^n\\ (3) 

n 
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which can be rewritten as 

P{f ^i)= Y.{^i\U'<{tf,ti)blal^\n){n\a,fbj,^U{tf,U)\^i). (4) 

n 

Baryon and lepton number conservations allow us to represent the summation over all electron and photon basis 
state compatible with the initial conditions by the insertion of identity, 

P{f ^i) = {^,\U\tf,U)bl^al^a,fhj,^U{tf,ti)\^,). (5) 

This expression is finally written as 

P(/ ^i) = Tv[OU{tf,ti)piU\tf,U)] (6) 
where the trace is taken over the full Fock space, 

Pi = &J,4l*9a«)(*sa«l«e,6p, (7) 

denotes the initial density matrix, and the Hermitean operator 

O = bl^al^aefbpf (8) 

handles the final state of the colliding particles. It is a well known procedure to sum transition probability over 
unresolved final states. But, the sum in Eq. (4) is over a too-large space and cannot be handled in the usual, explicit 
manner. The advantage of the form (6) is that it is an expectation value of a single observable incorporating the 
summation over the final states in Eq. (3). 



III. SYSTEM-ENVIRONMENT ENTANGLEMENT 



Before applying perturbation for the expectation value (6) , we point out the advantage of this formalism in finding 
the influence of an environment on the system studied. System stands for the collections of degrees of freedom with 
time evolution followed, environment denotes the remaining part of dynamics. The boundary conditions in time 
differ for these two components: the environment always follows open-ended, free time evolution encoded as an initial 
condition problem. The system may obey both intial and final conditions, like in scattering experiments. 

Let us introduce the generic field variables (/) and x for out system and its environment and the action S[(j), x\ 
governing the dynamics. The initial state is given by the density matrix Pi[{<p^ ,X^), {<P~ jX~)] where the fields with 
tilde depend on spatial coordinates only, (p = (f>{x), etc.. Our goal is to find the reduced density matrix 



(9) 



for the system. 

Consider first the complete density matrix subject of the time evolution 



Pt = U{t,ti)piU\t,ti). (10) 

The reduced density matrix is obtained by tracing out the environment suggesting the introduction of the generator 
functional 

where the trace is over the environment and T* denotes antitime ordering. The path integral representation of this 
functional, 



The boundary conditions in time, suppressed in this equation for the sake of better readability are the following. 
There is an integration over the initial field configurations with the weight factor Pi[{4''^ ,X^), {4>~ ,x~)\. At the final 
time we impose CTP boundary conditions for the environment, x+(t^,a;) = to assure an unconstrained 

time evolution for the system, the integration of the final configuration x corresponding to the trace operation in 
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the environment sector in Eqs. (9), (11). Finally, the system trajectories obey open time path (OTP) boundary 
conditions, (j)'^{tf,x) = ((y^{x). 

We write now the bare action of the theory as the sum of a system action and the rest, S[(f), x] = Ss[(j)] + Se[(l), x]) 
and consider the effective theory for the system by integrating over the environment variables, 

^iW\j;4>^] ^ j _fjj^jgiS3[<?i+]-iS.[</.-]-iS7[</.+ ,0-]+i3-<? 

where the influence functional Si[(j)^ [28] is deflned as 

The effective vertices, terms in the influence functional, can be classified as direct and entangled, the latter coupling 
<j)+ and <j)-, Si[<j)+,<j)-] = - + 5e [</>+, </>"], Se/Sip+Sip- ^ 0. In the OTP scheme interactions are 

represented by couplings among fields on the same time axis, eg. the interactions within the system give rise 5'd[(/)^]. 
Each contribution of the path integral (12) for a given pair of environment configurations x^ gives a factorisable 
contribution to the density matrix, (i)^; x"*", X ] which represents a pure system state. The CTP boundary 

conditions for the environment couple x^ with x^ and generate 5*6 [(/>^, (p~\- It is the integration over the environment 
configurations in Eq. (14) that produces a density matrix of a mixed state. 

The path integral representation of the reduced density matrix makes two phenomena related to the classical 
limit, consistency and dccohercncc, particularly explicit. The smallncss of the off-diagonal matrix elements, the 
linear superposition of different coordinate states, is a measure of decoherence [29], the quality of ^ as a pointer 
variable [30]. Decoherence is induced by the interactions with the environment. Let us now consider a pair of system 
trajectories 0^, their contribution to the reduced density matrix is the integrands on the right-hand side of Eq. (13), 
which is suppressed due to interactions with the environment when Sj5'/[0+,0~] <C (fi- = 1). Hence, the more 
negative is the imaginary part of the influence functional, the more decohered is the pair . 

Decoherence characterizes the reduced density matrix at a given instant of time. It builds up in time and this 
dynamical process is called consistency [2-4]. Let us consider now two trajectories </>"'" and 0~, <p^{tf,x) = (j)~{tf,x) = 
^{x) contributing to the diagonal matrix element in Eq. (13) at (p. There are four contributions, four pairs or 
trajectories related to 0+ and (p~ , namely (p = (</)+,(/)+), ((/>+, i/i"), ((/>",(/)+) and {(p~,(f)~). This pair of trajectory is 
called consistent if their contribution to the probability of finding (p in the final state is additive. This situation is 
approached when the contribution of (p = {<p'^,<p~) is small, QSilcp'^ ,(p~] <C 0. 

Relation between decoherence and consistency is clear: A pair of trajectories is consistent if the imaginary part of the 
infiuence functional of the CTP formalism receives sufficiently negative contributions from the time t' of propagation 
when these trajectories are different, (p'^{t',x) ^ (p~{t',x). But stopping the time evolution at such a time and 
considering the contribution in OTP scheme this integrand represents a contribution to the reduced density matrix 
Ps[4<^{t' ,x),(p~{t' ,x)] and indicates decoherence at that particular time. In other words, decoherence sustained in 
sufficiently long time makes consistency. 

IV. TRANSITION PROBABILITY 

We introduce now the reduction formula for the electron-proton elastic scattering, together with its perturbative 
evaluation. 

A. Reduction formulsis 

The transformation of the transition probability of an inclusive scattering process, (4) into an expectation value 
(6) requires the use of the reduction formulas [31]. It is obvious to generalize this scheme for the CTP formalism and 
the result can be easily presented in terms of the generator functional for connected Green functions, constructed 
by introducing external sources coupled linearly to the elementary fields in the Lagrangian. To handle both time 
evolution operators in Eq. (6) one introduces independent external sources for U and indexed with + and -, 

respectively, r)e = {ilt^V7) ^^'^ Vp — iVp^Vp) for the electron and proton fields and j — for the gauge field. 

The generator functional is, therefore, defined by 

^w\j,fi.,fiA = TT[U{tf,U;r,+ ,fj+,j+)p,U\tf,U;-7j.,-fj.,-r)] (15) 



(13) 



(14) 
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{h = c=l) where 

U{tf ti-rj fj j) = 7"[e~'-/t*/ ■'^''^^[^^^^"^^(''•^(^^'^•^(^'+'^^(^''''^^^^^"^''(^^'^''^^^1] (16) 

denotes the time evolution operator in the presence of the external sources rj,f],j, and H(x) stands for the energy 
density, r = e or p. Each degree of freedom appears twice, once in the time ordered exponentialized form oi U(ti,tf) 
and another one in its Hermitian conjugate, involving antitime ordering T*. The path integral representation of (15) 
is 

where the integration is over pairs of trajectories ipr = {ip^,tp~), A = {A^,A~) and the action 

S[A, i,,,i,,]=Y,j dxMx){G-^ - erKx))Mx) + \ j dxA{x)b^^ A{x) (18) 

contains 4= A^^i^ with the Dirac doublet matrices 7** = (t**, — (l**)*)- The fermion inverse propagators are 

A-i _ (ip-rrir + ie ^ , 



-70(i^^ + ie)t7o; ^""scr 

nie = m and rup = M being the electron and proton mass, respectively. The inverse photon propagator 

^0-;. = 9,. (° -□ + ie) + ^b'c, (20) 
is given in Feynman gauge. The boundary conditions in time are 

Ki^f,^) = A-{tf,x), (21) 

due to the trace in Eq. (15). The closing of the path at the final time can be implemented on the level of the free 
propagators. Since the free propagators can be most easily derived in the operator formalism the actual form of G^q 

and Dgc ^^^^ ^'^^ important for us. 

One has to introduce a gauge invariant cutoff, for instance, dimensional regularization after Wick rotation, and the 
corresponding counterterms in the exponent of Eq. (17). Since the boundary conditions in time do not affect, the 
UV divergences the renormalization and the counterterm structure are identical with the single time axis formalism 
and will i)e suppressed in what follows. To remove IR divergent tadpole contributions we assume the presence of a 
classical, homogeneous external charge density that neutralizes the electron gas. 

We can finally turn to the transition probability (6), which can be obtained by means of the generator functional 
introduced above; by repeating the steps followed in deriving the reduction formulas. The transition probability of 
the elastic scattering process e{pi) + p{p2) — >■ e(gi) +^(92) reads as 

P{f^i) = Z-'^Zp'^ j dxtdxtdytdytdx^dx^dy^dy^e''i'^yi-y^^+''i^^yi-y^^-'P'^''i-''^^^^^^ 



(22) 



where Z denotes the wave function renormalization constant, expressing the proportionality of asymptotic and in- 
terpolating fields, Ve = Zl^'^t/j^^ and tpp = Zp^^ip^^. The Dirac spinors uj and uj describe the initial and final 
states, respectively, for electrons {j = 1) and protons {j = 2). This expression will be evaluated below by a partial 
resummation of the perturbation series. 
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The calculation of scattering probability brings out an important difference between the TFT and the CTP formal- 
ism. In TFT, the forward-pointing time axis is used to follow physical processes and the other time axis that contains 
the so-called ghost degrees of freedom is a formal device to calculate the thermal averages with real time dependence. 
The corresponding generator functional for QED is 



ZTFT = Tr[U{tf,ti;rj+,f}+,j+)pl-'U^itf,U;-rj-,-fi-,-j-)pt], 



(23) 



where pi is the density matrix in thermal equilibrium and < s < 1. It is easy to see that the expectation values 
obtained by functional derivatives with respect to the external sources 77+, 77+ and j"*" evaluated at f) = ^ = j = are 
independent of s, whose value can be freely set to simplify the free Green functions. One can obtain in this manner 
the thermal average of an observable or the transition probability of a process as long as equilibrium prevails, namely 
the condition [pf,U{t,t';0)] = holds after carrying out the observation or completing the process. But, it may 
happen that the measurement or the process generates nonequilibrium phenomenon such as the decay of the "vacuum 
polarization cloud" in the gas, made by the colliding particles. In such a case, we need s = 0, the choice that makes 
the functionals (15) and (23) identical, renders the ghost degrees of freedom physical and brings TFT equivalent with 
CTP. Another special feature of CTP, the choice s = 0, is that it makes the reduced density matrix available and 
opens the way to address the issue of quantum-classical transition in expectation values or in processes. 



B. Perturbation expansion 

The free generating functional, 
Wo\j 

contains the CTP propagators with Fourier transforms 

Go. (5) = j dxGor{x,0)e'''\ 
Mq) = J dxDoix, 0)6^1=', 

given by 



(24) 



27ri5(g2 - TO2)e(-gO)^ 



2m5{q^ - ml)@{q°) ^ 



-|- 2'Ki6{q^ — m^)n, 



q''—7n^—ie / 

where the occupation number at zero temperature and finite density is given by 

rigr = e{q°)e{pr - e^r) + e{-q°)Q{-egr - Mr) 



1 -1 

<ir[_i 1 



(25) 



(26) 



(27) 



with /ie = /U, /Up = and eqr = \/ -|- [32]. The photon propagator is used in Feynman gauge, = g'^'^Do{q) 

with 



27ri(5(92)e(-gO) 



The structure of these matrices reflects a general rule, valid for any local composite operators. 



iD{x,y) 



{T[^{x)c^{y)]) {cb{y)cl>{x)) 

{mm) {T[mm]Y 



for bosonic operators and 



iG{x,y) 



[ -(0|V(a;)^(y)|0) (0\T[{j°^{yMmh')]\Oy 



t[Px + rn) \ _Qf _Gnj^iQ' 



(28) 



(29) 



(30) 



for fcrmionic operators. Furthermore, D^{q) = D"{—q), D^{q) = D^{—q), G"{q) = G"{—q), G^{q) = G^{—q) are real 
and {q) = —D^{—q) and G^ [q) = —G^{—q) are imaginary. Note that the off-shell contributions are in the real 



7 



part of the diagonal blocks only. The environment contributions are always on shell. In case the photon propagator 
D" and are the near and far field, 



give the retarded and advanced Green functions of classical electrodynamics. 
We return now to the generator functional (17) which assumes the form 



(31) 



o-=±l 



dx- 



iWo\j^'>le,fie,flp,flp] 



The lowest order contribution to the transition probability (22) comes from the 0{fjprip f]'^r]ge^) term 



(32) 



4! 



a dx 



a=±l 



4! 



Sj'^t' {x) 
dxdyne{x)Goe{x, y)fie{y) 



' H) 



222! 



dxdyj'^{x)Don„{x, y)f{y) 



4! 



/ 



dxdyr]p{x)Gop{x, y)fip{y) 



(33) 



in Eq. (32) where we can use Z,, = Zp = 1. 

Our system consists of the colliding electron-proton pair and photons, their environment being the remaining 

electrons in the gas with occupation number (27). The operators of the free equation of motion that truncate the 
Green function in Eq. (22) when they act on the external leg according to the reduction formulas extract the residuum 
at the mass shell. Owing to the identity 



5{x)x = 



(34) 



this operation suppresses all mass-shell contribution, such as the off-diagonal CTP blocks and the 0{nk) environment 
piece of the electron propagator (26). In physical terms, we find the usual OTP propagator handling our system, whose 
asymptotic sector contains the colliding particles only. Naturally, the soft photon component of asymptotic charged 
states [33], ignored in this work, should be constructed as in the usual scattering processes without environment. 

There are two topologically different graphs characterizing the contributions in (33), they are depicted in Fig. 1. 
Graph (a), 

-e^J durdu2D++{u, - U2)G++,^{y+ - «i)7,^,,+ G++^+(«i - xt)G++^,^iy+ - uM+,,+ G+t,An2 - 4) 

dvidv2D^^^,{vi - V2)G-r^,_{y^ - vi)-('^,l^,_G--^_{yi - x^)G-I^,_{y^ - V2)'yd'-b'-Gb'-b-i^-^ ' a^^)-(35) 

contains the CTP-diagonal contributions of the photon propagator and displays interaction betwcxai the colliding 
charges by exchanging a photon. The state of the colliding particles remains pure after interacting and the transition 
probability is the amplitude times its complex conjugate. Graph (b) in Fig. 1, 

-e'l du,dU2D+-^,{u, - V2)G++^,Ayt - «l)7e^„'+G++„+(ni - 4)G+V(y2+ - «2)7d'+6'+G'++ +(^2 - 4) 

X J dv^dv2D-+Jv, - U2)G--,_(yr - viH:-\,-G--,-{vi x-)G-r,,_{y^ - V2)Y,.\,-G--^_{v2 - x^), (36) 

includes the CTP off-diagonal part of the photon propagator. Note that these blocks contain the one-shell Wightman 
function and have to be represented by a line of the corresponding particle, namely, the photon that connects the two 
time axes by passing between them at the final time. The graphs, therefore, are equipped with periodic boundary 
conditions in time; lines reaching the final time at one time axis enter at the other time axis on the other side of the 
diagram. Therefore, this graph describes the contributions of final states containing two photons to the scattering. 
When a photon is attached to the same particle at a different time axis, then this contribution represents the soft 
photon component of the colliding charges. But, we ignore forward scattering and the leading-order graph must 
contain photons that are exchanged between the colliding particles. Since the support of the off-diagonal CTP blocks 
of the free photon propagator is on the mass shell, the emission of a single photon is kinematically excluded and this 
unusual graph is vanishing in the absence of the electron gas. 
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FIG. 1: The leading order graphs contributing to the transition probability: (a): direct, (b): entanglement contributions. The 
circle represents the initial density matrix. The CTP-diagonal block of the photon propagator connects charged particle lines 
on the same time axis in graph (a); the off-diagonal blocks couple the two time axes and their photon lines traverse the final 
time, the left and right end of graph (b), with periodic boundary conditions in time, i.e., lines ending at the same altitude at 
the left and right end of the graph are joined. 



t-channel resummation 



The higher-order corrections to the transition probabihty (22) can be partially resummed by means of the skeleton 
expansion, where graphs with diflFerent self-energy or vertex insertions are resummed. The result is a graph similar 
to the lowest-order contribution except that the propagators and vertices arc replaced by their exact expressions. 
Such a resummation represents an important improvement for kinematical regimes where one of these quantities 
assumes values strongly different than in leading order. We shall be interested in scattering process where collective 
particlc-hole excitations take place in the t channel. For this end, the usual resummation of the photon self-energy is 
necessary. 

A natural consistency requirement for any partial resummation is to preserve gauge invariance to properly separate 
the physical sector from the IR singular, nonphysical gauge modes. A simple way of checking gauge invariance of a 
given calculation of the irreducible vertex functions is to construct the effective action and verify its invariance under 
gauge transformation. Photons arc electrically neutral therefore their effective CTP action obtained by the partial 
resummation of the self-energy S is the free Maxwell action 



except that the inverse free photon propagator Dq^ is replaced by the dressed one, 

-E. 

The self-energy is transverse, 9^S^y(x,2/) = and our partial resummation is consistent. 
The reduction formula for the transition probability (22) yields the sum 

P{f ^i) = Pd + Pe, 



(37) 



(38) 



(39) 



where the first and the second term corresponds to expressions (35) and (36), respectively. The resummation of the 
phonon-self-energy amounts to the replacement of the free photon propagator with (38) in Eqs. (35) and (36). The 
CTP-diagonal direct term. 



Pd = [ dxldx'^^dyldy^e 



+ 

+ 6+ 



(m2 - X-^) 



X j du,du2D++{m ~ u,)G++^,^{y+ - u^h^,^^,^G+t^Ani - xt)G++^,^{y+ - u,)r^,+,,+ G+ 
X j dvidv2D-,-,{vi - V2)G;--,_{y^ - vi)j^!\,_G--^_{vi - x^G".- ,_ (2/3" - V2)-fa,-f,,-G--^_{v2 - x^) 



(40) 



contains the contributions of the colliding particle interaction. Since the support of the off-diagonal CTP blocks 

is the mass-shell all free photon lines in the Schwinger-Dyson resummed graphs stand for D^^. Hence, the number of 
2=*=^ insertions is always even in P^. An 0(i;++(I]"' )^) piece of the geometrical series, resummed Eq. (38), is shown 
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FIG. 2: Two higher order graphs contributing to the transition probabihty: (a): direct, (b): entanglement contributions. The 
particle-hole pairs dress the propagation of the photon in each time axis separately in the direct graph and traverse the final 
time in the entangled graph. 



in Fig. 2 (a). The self-energy insertion S++ describes a particle-hole pair creation contributing to the interactions of 
the colliding particles on the -|- axis. The correlation of the two colliding particles on the — axis represents interaction. 
In fact, each particle interacts with a particle-hole pair that does not interact with each other but forms a final state 

that is matched with two particle-holc pairs of the other time axis. These pairs arc the source of correlation because 
they interact by arising from the same photon. Therefore, the matching of the asymptotic states at the final time 
can transfer correlations from one time axis to the other one. Note that the off-diagonal CTP blocks of the photon 
self-energy, E^^, have a special role in the transition probability; they represent the asymptotic particle- hole state 
contributions, the backreaction of the colliding charge on the electron gas. For instance, the n-fold insertion of the 
one-loop Ti^^ captures the dynamics of n asymptotic particle-hole pairs. 
The contributions of the CTP off-diagonal blocks of the photon propagator, 



(T = ± • 

xuifiipy+ - m)u2f{ipy+ - M)ulf{-ipy- - m)ulf{-ipy- - M) 



X j duidu2D+-,{ui - V2)G++^,+ {y+ - ui)'y^,+^,+ G++^+ (wi - x+)G++^,+ {y+ - «2)7d'+6'+GjV («2 - 4) 

X j dvidv2D-+{vi - U2)G-r^,_{y^ - vi)-i''^,l^,_G--^^{vi - x^)G-^^^^ {V2-X2) 

A-iP^t - - M)u2^{lp^-m)uU^p^- - M)ul,, (41) 

reflect colliding charge-gas entanglement. An 0((5]+~)^) graph contributing to is depicted in Fig. 2 (b). The 
colliding charges produce particle- hole pairs on both time axes. The matching of these pairs at the final time establishes 
a coupling between the time axes which qualifies entanglement according to the remarks made after Eq. (14) . 
A straightforward calculation gives the transition probabilities 

= yW(27r)4%i-pi+g2-p2)|rdP, 
Pe = V^^\2TT)H{q,-Vi+q2-V2)\Te\\ (42) 

where V^^^ is the four- volume and the "transition amplitude" squares are given by positive semidefinite expressions 

\Ta? = e4|(,liy.7^,i,)(,l2y7''«20^J.+ ('Zi-Pi)l', 

\%? = e4|(«l/7^«l^)(«2/7'^"2i)*i)+-(gl-Pl)|^ (43) 

by means of the relations D'j^'J (p) = [D~^'~'^ {~P)]* for the CTP propagator for any, not necessarily free, local 
four- vector operator. 



D. Photon self-energy 



Before presenting the photon self-energy let us work out two useful parameterizations of Lorentz tensors appearing 
in this work. The photon self-energy S^'' and propagator D'^'^ have two intrinsic Lorentz vectors, the four-momentum 
and a unit timelike vector, assumed to be = (1,0) in the rest frame of the gas. The symmetric tensor M, 
denoting either self-energy or propagator, constructed by the help of these vectors has three free scalar parameters, 
for instance the coefficients of the terms q'^q", n^n^ , and q^n" + n^q". Current conservation allows us to ignore the 
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longitudinal part of the photon propagator in Eqs. (43) hence transversality, q^M^^{q) = 0, reduces the number 
of free scalars to two that will be chosen Mg = M'^'^g^f^ and M„ = Ufj^M^'^Uv This parameterization is useful to 
calculate loop integrals, cf. Eq. (47) below. 

When the inverse of a tensor is needed for the calculation of the propagator by means of the Schwinger-Dyson 
equation (38), then the parameters Mg and M„ are not useful anymore. Instead, it is more advantageous to split 
the four-dimensional transverse subspacc, identified by the projector T*"^ = gi^^ — q^<f jq^ into three-dimensional 
transverse and longitudinal subspaces, T = + Pt, where the projectors 

p^. _ _ /O 0\ 

= {1 7i) («) 

corresponding to the transverse and longitudinal modes, respectively, in Lorentz gauge are given in terms of the 
three-dimensional longitudinal and transverse projectors, L = n®n,T = \ — L, with n = q/|q| and v = q^/\q\. 
The inverse of the matrix 



M^'''{Mf,, Mt) = MiP^" + MtPt 



(45) 



can be obtained in an obvious manner within both subspaces, M ^{M(, Mt) = M^^'^{l/M(, 1/Mt). It is easy to check 
the relation 



Me = {l- z/2)M„, Mt = \[Mg + {v^ - 1)M„] 



(46) 



between the two different parametrization. 

The one-loop photon self-energy in an environment [34] has already been studied thoroughly; we need its CTP form 
[35] for a zero temperature nondegenerate electron gas of Fermi momentum k, 



(47) 



The 2x2 CTP blocks are given by the matrices which can be written as sums of the vacuum and electron gas 
contributions as [32] 



— vac ~1~ ^x gas 

where x stands for the letter g or u and the self-energy displays the structure 

[n^ + iw -n" + iir 

with 



1 

-1) ■ 



The vacuum contributions contain the functions 



KaM = 



37r 



2 1 



2m' 



I Am? 



— 1 arccot< 



-1-1 



ia\ m 

2 



nL(9) = 

with a = e^/47r and uj{q) = \/rr? +p^, 



u/{p) 



3|9| 



n^.„c(9) = 35in^„,(9), 

Kn" (a) 

u vac\^) ^2 ' 



(48) 



(49) 



(50) 



(51) 



(52) 
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and similar expression for 11^ We shall consider below the nonrelativistic regime, \q^\ < m where n^^p(g) = 
and ^K^M = 0- 

A zero temperature nonrelativistic gas of electrons with Fermi momentum k <^ m gives 

Kgasil) = nS,„,(g) + nS,„,(-g), 



with 



and 



X gas 
^■^x gas 



(9) 
(9) 



K gas{<l)+K gasi-Q), 



n" 

g gas 



(9) = 



11" 



M) = 



H gasiQ) 



HI 



2ak m 
7r|qr| 

2ak'^m 
7r|q| 

ak'^m 

ak^m 
ak^m 



Hi). 



n 



u gas 



(9) 



2 777, 2 
r,2 



M{q), 



q 



2m? 

„2 



q 



Am? 



Niq), 



- m, 

m ' 



L{q) = r+-(l-r2)ln 



M{q) = e(l-|r|)(l-r2), 



N{q) 




\q\>2k, 
\q\<2k, 



1 < r < 1 



kl 



< r < 1 



l9l<2fc, -i<r-<l-M 
The q dependence of the functions L{q), M{q) and N((i) is given through the combination r = {(f 



(53) 



(54) 



(55) 
2mq^)l2\q\k. 



V. PHOTON PROPAGATOR 

The actual resummation of the geometrical series with the self-energy is carried out by the algebraic solution of the 
Schwinger-Dyson-Kadanoff-Baym equation, Eq. (38). We first consider the CTP and after that the Lorentz tensor 
structure of the propagator. 



CTP blocks 



The CTP blocks structure of the propagator is shown in Eq. (29) and its inverse, both D^"^ and H, have a common, 
slightly different 2x2 block structure of Eq. (49). related by the change of basis by a. The Green- function-like 
structure is inherited during multiplications where a is treated a metric tensor, i.e., the Green function and its inverse 
are interpreted as covariant and contravariant tensors, respectively. In fact, one can easily prove the relations 



{DiaD'2a---Djf 
{DiaD2a ■ ■ ■ DjY 



D^D^---D'^, 

{D^aD^a-- ■ Dj.^yD] + {D^aD^a ■ ■ ■ D^.^fD), (56) 

for a set of Green function like 2x2 matrices, Di,. . ., Dn by recursion [32] . These relations are useful when the 
propagator (38) is written as 

1 



D = 



1 - Doana 



-Do 



(57) 
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and the first equation in (56) gives 



= — D^^ = -. (58) 

1 - £>gn« £>o ^ - n° 



The imaginary part is obtained by using the second equation in (56) 

1 



1 - Doalla 



i=i 3=0 

^ ^ ^gn°+ , l^^A DoaUay, (59) 



1 - Doalla J " 1 - D^^Ur 



yielding 



1 - Doalla 



= (DiW+DT.W) (60) 



and 

D' = D'^WD" (61) 

where the equation DqDq^ = 0, an apphcation of Eq. (34), has been used. Finally, Eqs. (58) and (61) together with 
(29) and (31) give the CTP structure of the resummed photon propagator. 

B. Lorentz blocks 

Next, we work out the retarded and advanced CTP blocks of the Lorentz tensors appearing on the right-hand side 
of the Schwinger-Dyson equation, (38). The inverse free propagator is usually obtained either by the identification 
of the kernel of the quadratic part of the action or simply by the inversion of the propagator in momcntiim space. 
Neither of these ways is trivial in the CTP scheme. The two time axes are coupled at the final time in the path 
integral (17). One cannot even have an inverse propagator that is invariant under translation in time. The CTP 
propagators can nevertheless be derived in a straightforward but rather lengthy way in the path integral formalism by 
carefully implementing the boundary conditions in time for free fields at finite At and performing the limit At — )• 0, 
followed hy t f — ti ^ oo . A translation- invariant kernel and diagonal form in frequency space is recovered without 
problem after the second limiting procedure. 

A shorter and more promising way is to derive the propagators (26) and (28) in the operator formalism where they 
are automatically diagonal in the Fourier space in the limit tf — ti —>■ oo. But, the difficulties are encountered in this 
case when the inverse of distributions, handling the mass-shell singularities, is sought. A simple, natural way out is 
to regulate the distributions by the replacement 

= (62) 

and performing the inversion with small but finite e. The result one finds in this manner is 

Note the transmutation of the coupling between the two time axes. One the one hand, this coupling is e- independent, 
0{e°) and is localized in time at t = tf in the action written in the space-time for finite tf — tf, like in Eq. (18). 
On the other hand, when we rewrite the action in the momentum space after the limit tf — ti oo by means of 
the free inverse propagator (63) then the boundary conditions in time arc difficult to trace and the coupling between 
the time axes is represented by a weak, 0{e) term in the free action which acts during the whole time evolution in a 
time-independent manner. 

The retarded and advanced part of the self-energy (47) are 

nf = (l-^.2)^i uf = l[nl + i,.^-i)ul] (64) 



13 



according to (46) where 11° = 11" ± IT-''. An important feature of the self-energy component 11* is that it is always 
nonpositive and we can in this manner ignore the second, 0{e) term of the free propagator (63) beside the self-energy, 
leading to Dg"^" = -q'^{Pe + Pt) and D'^^ = D^^' = 0. 

The Schwinger-Dyson resummed retarded and advanced propagators are found first by means of Eq. (58), 

Dt = ^- Dt = , . (65) 

(t.2 - i)(q2 _ ns) in« + (z/2-i)(ins + fe2) 

The imaginary part of the propagator, given by Eq. (61) turns out to be 

Di = DlD^Ul,{l-,.% Dl = ^DlD1[Ul + {i^'-l)K]- (66) 
This propagator inserted in Eqs. (42)- (43) gives the partially resummed expression for the transition probability. 



C. Nonrelativistic spectral representation 

We close this section by a brief summary of the spectral representation without relativistic symmetries for the two- 
point functions of local transverse vector operators, in particular for the photon field or the electric current within the 
CTP scheme. The spectral representation of a relativistic propagator is an integral over a Lorentz scalar where the 
integrand is the product of a spectral weight and a free propagator containing this Lorentz scalar. One expects similar 
integral representation in the absence of Lorentz symmetry as well, except that the integral variable would not be 
scalar. An important advantage of the spectral representation is that it provides a common parametrization of causal, 
retarded, and advanced propagators. This feature will arise in a specially clear manner within the CTP scheme, the 
natural basis to introduce these different Green functions. Another issue to detail is the way the parametrization (45) 
can be carried over the spectral functions. 

Let us suppose first that we have a local real field (f){x) with propagator given by Eq. (29). The spectral function 

27rp(^)(g) = J2i'^n)'^i<l-Pn)mm\n)\' (67) 

n 

is a sum over the eigenvectors |n) of the Hamiltonian with four- momentum p^, = 0, and > for n > 0. The 
spectral function is vanishing for negative frequency. 
We introduce the CTP matrix 



Vo[iv,oj ) - -c/ ,,2 _ „,/2^c,^.,^ 1 (6«) 



uj^—oj'^—ze 



and elementary steps followed in the derivation of the usual spectral representation yield 

b{q)= / du:'^Do{q\J)p{u:',q). (69) 







In case of rotational invariance, the spectral strength is a two-variable function, p^'^\q) = p*^ ■*((?" , |<7|). Note that the 
spectral integral (68) is trivial for the CTP blocks and the spectral function can be read off directly from the 

propagator since 

2TTp^'''\q) =iD-+{q) (70) 

holds for q° > 0. In case of a free relativistic field of mass m we find p'^'f'^q) = 5{q^-m^). 
The generalization of the spectral function for a local vector field An{x) is a Lorentz tensor 

2^pi^\q)^Y.^2^)'^^'l-PnmMnn){n\A,{m (71) 

n 

and will be characterized by the parameters 

p^^^ = {l-v^)u<'u''p^,, 

P^^^ = \[9'"' + {^^-^)u^u'']p,, (72) 
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in case of the gauge field. This choice leads to a useful factorization of the CTP and the Lorentz indices, 

D^'^iq) = De{q)Pr + Dt{q)Pr (73) 

with 

De{q) = / du'^bo{q^J)p[''\uj',q), 
Jo 

POO 

Dt{q) = / ckv'^Do{q'>,u;')pi^\u>',q). (74) 
Jo 

The spectral functions can again be obtained from the Wightman function, 

iD-+ = 2'Kp[^'^Pt + 27rp^^^Pt (75) 

for > 0. 

VI. ELECTRON-PROTON COLLISION 

The calculation of the collision process in the double time scheme gives not only expectation values and transition 
probabilities, but provides information about consistency and irreversibility, necessary conditions for classical limit. 

A. Cross section 

Once the transition probability, displayed by Eqs. (42)-(43) for the scattering e(pi) + p{p2) — > e.{qi) +^(92), is 
found, then the construction of the cross section is straightforward and standard. Averaging over the initial and sum 
of the final fermion polarizations in Eqs. (42)- (43) yields 

\Td? = ^^^{2Sft[(piD++(r)p2)(gii?++*(r)?2)] +25i[(pii)++(r)92)(?ii5++*(r)p2)] 
+tr[£)++(r)£)++*(r)](m2 -pigi)(M2 -^292) 

+2{M^-p2q2)^\piD++ir)D++*'^{r)qi]+2{m^-p^q^)^\p2D++*^r)D++*ir)q2]}, 

\%? = ^^^{2SR[(piI)+-(r)p2)(9il)+-*(r)92)] +25R[(pi£>+-(r)g2)(gii?+-*(r)p2)] 

+ir[D+- {r)D+-\r)]{m' ~ Piqi)iM^ ~ P2q2) 

+2{M^ - p2q2W\piD+-ir)D+-^r)qi] + 2{n? - pigi)K[92£'-+^(-r)r'-+(-r)p2]}, (76) 
where r = q\— p\. The cross section is finally obtained in the usual fashion, 

I 



where the intensity of the transition, 



(77) 



, d^qid^q2mM 4 2 /-7Q^ 

I (2^0M^IM^(2") S{q,-p,+q2-p2)\T\ (78) 



with w(g) given after Eq. (52) and ^{q) = y/M^ + q^. The integration is over the region S in the final momentum 
space covered by the detector. 

The cross section (77) is manifestly Lorentz invariant in the vacuum, fc = 0. In the presence of environment, the 
electron gas, it depends on the average momentum of the environment. When the gas is at rest in the center-of- 
mass frame of the colliding particles, then the exchanged photon has vanishing energy. This makes the entanglement 
contributions to the cross scattering, |7^P, vanishing. 

To see the effects of entanglement, we choose both the colliding electron and the electron gas at rest in the laboratory 
frame. In other words, we consider an experiment where an electron gas is taken as target into a proton beam and 
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(a) (b) (c) 

FIG. 3: Logarithm of the difTerential cross sections (a) log(T„, (b) logaj and (c) log(atj + a^) as functions of P/k and Q/tt. 

those scattering events are filtered out where the initial electron is in the state p — 0. The differential cross section 
in this frame, corresponding to a given beam particle momentum p2 and scattering angle dfl, is 

^^l^^l -irr (79) 



d\p2\dn 167r2|p2|(m + n{p2) - cj(qi)) ' 

Cross sections, calculated by using the transition probability |7dP or |7^P in this expression will be denoted by ad 
and CTe, respectively. 

The natural parameters of the cross section (79) are the beam momentum and the scattering angle. But the 
kinematical regime allowed on this two-dimensional plane is not a simple rectangle. To simplify the independent 
variables, we bring the system into the center-of-mass frame for the colliding charges by performing an appropriate 
Lorentz transformation. Since the coefficient of the transition probability on the right-hand side of Eq. (79) is Lorentz 
invariant, all we need to do is to parametrize the scattering amplitude by the natural variables of the center-of-mass 
frame, namely, the charged particle momentum P and scattering angle Q. The particle momenta in the center-of-mass 
system, written as 

Pi = (r2(P),-P,0,0), Qi = (f7(P),-Pcose,-Psine,0), 

P2 = (cj(P),P,0,0), Q2 = {uj{P),Pcose,Psme,0), (80) 
is used to express momenta in the rest frame of the electrons, 

P^ - vP" 

P 



where the boost velocity in the x-direction is w = P/uj{P). We arrive in this manner at the expressions 

P + V^p 



Pi 



/Pcos6 — vflp ^ , , , 

q2 = [—^===^,Psme]. (82) 



The cross section (79) now reads as 

da^ _ \qi\{n{P)+oj{P)) 



dPdVL 167r2 \p2 \{m + n{p2) - uj{qi)) 



p2 



f^(P)w(P) 



\r\'. (83) 



B. Diagonal contribution 



The numerical results, presented below have been obtained for k — 0.003m, corresponding to 2.25eV Fermi energy, 
typical metallic density at the lower edge of applicability of perturbation expansion. 
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FIG. 4: Spectral strengths (a) 2-Kpf of the electric current, (b) 2-Kp\ and (c) 2iTp\ of the photon propagator as functions 
of \q\/k and z = mq° /k^ . The additional curve belongs to = O.Stt. 




5 2.0 



FIG. 5: (a) _D" and contourplots of (b) D", (c): 2Tip\ as functions of \q\/k and z = mq /k . The line Q — 0.5n is shown, as 
well. 






FIG. 6: (a) |Td|^ (b) 27rp<'*', (c) H?, (d) 2-Kpf >, (e) \Td\'' + \Te\'' (higher curves), ITdl" (lower curves), and (f) \Te\'' as functions 
of P at O = O.Stt (left curve) and 9 — 0.5n (right curve). 
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The differential cross section ay in the vacuum, fc = 0, is shown in Fig. 3 (a) as the function of dimensionless 
center-of-mass momentum P/k and scattering angle Q/n. It diverges when P — >■ or — > 0, but its value is cut off 
in the plot for clarity. The cross section, ad, depicted in Fig. 3 (b), shows two obvious effects. One is screening, a 
suppression for small energy exchanges. The other is a valley and a rim appearing at a certain exchange energy. 

It is instructive at this point to look into the spectral strengths. The spectral strength 2ttp^/\ introduced in Eq. (72) 
is a measure of the amount of particle-hole states available for the longitudinal electric current operator. It is shown in 
Fig. 4 (a) as the function of the dimensionless momentum |<7|/fc and frequency z = mq^ /k"^ . It is nonvanishing around 
the free-particle dispersion relation and displays a steep, linear increase at low frequencies within the momentum 
range \vq\ < 2k ending at a sharp edge. This and the transverse spectral strength of the current arc built into the 
photon propagator in a nonlinear manner, as indicated by Eq. (65). The resulting spectral strengths for the photon 
propagator are given in Figs. 4 (b) and (c). 

Another curve is shown in the figures, as well; it is the projection of the line parametrized by the momentum P in 
the center-of-mass system at 6 = O.Stt. The nonlinearity of the Schwinger-Dyson resummation transforms and shifts 
the maximum of Fig. 4 (a) towards larger frequencies in the photon spectral strength seen in Fig. 4 (b). The curve 
= O.Stt reaches its maximum at \q\/k ^ 1, at the local maximum of ad in Fig. 3 (b) at the same scattering angle. 
The spectral weight of transverse photons, shown in Fig. 4 (c), displays a narrow diverging peak at the origin and the 
line O = O.Stt slides down on it as P increases. This singularity is suppressed by the momentum variables multiplying 
the propagator and gives small contribution around the rim in Fig. 3 (b), which is due to the increased number of 
particle-hole states to mix with longitudinal photons. 

The electric current spectral strength display maximum around quasiparticles, collective particle-hole excitations, 
defined by the root of ^{D^~^)~^ on the energy-momentum plane. The function D"^ and its contour plot are shown 
in Figs. 5 (a) and (b) together with the line 6 = O.Stt. The propagator is divergent along the heavy line starting 
approximately horizontally at z ~ 1, corresponding to plasmonlike collective excitations [32]. The imaginary part 
of the one-loop propagator is vanishing in this region; these collective modes decay slowly into several particle-hole 
pairs. $Sn^ starts to become nonvanishing when we arrive at the end of the solid line at 2; ~ 1.5. The lowest line 
separating different colors starting from this point on the contour plot runs approximately along the line ^{D^~^)~^ = 
corresponding to the zero-sound excitations. The longitudinal propagator, D", is finite along this quasi particle line 
where 911^ 7^ making the zero-sound strongly damped by single particle-hole pairs. The peak of the cross section 
at P ^ 1 appears at \q\/k ^ 1.4, approximately the closest point of the end of the singular solid line to the curve 

Q = O.Stt. The contour plot of pf^\ shown in Fig. 5 (c), indicates that the number of available longitudinal photon 
states reaches a maximum at the end of the solid line of Fig. 5 (a), where the plasmon and zero-sound lines meet. 
The lesson of these results is that quasiparticles around the Fermi surface, \q\ fc, 2; 1, are responsible for the rim 
in the direct cross section ad- 

To better locate the rim and valley of ad in Fig. 3 (b), two lines, one with 6 = 0.3tt and another with 6 = O.Stt, are 
followed in Fig. 6 as functions of the center-of-mass momentum P in different terms of the transition probability. It 
is better to look into the transition probability rather than the cross section because the former has no kinematical 
factors besides the expectation value (3). The transition probabilities at 6 — O.Stt and 9 = O.Stt are shown in Fig. 6 
(a). Their peak is approximately at the maximum of p^/^\ shown in Fig. 6 (b), supporting the remark made before 
about the coincidence of the maximum of the cross section and of the number of longitudinal photon states. The valley 

of the cross section seem to agree with the position of the maximum of p^^\ shown in Fig. 6 (c). This is rather natural; 

■ (j) 

faster-decaying particle- hole states make weaker transition probability. The sharp rim of pY separates two different 
kinematical regions, the last two cases in the definition of N in (55). This makes a characteristic, continuous but 
nondifferentiable singularity for the transition probability. It is the same location where the self-energy 11", depicted 
in Fig. 6 (d), has inflection point. 

C. Entanglement 

We now inspect the entanglement contribution to the cross section. According to the remark about the parity of the 
number of T,^^ factors in ad and ag whenever the latter dominates the former the colliding particle-gas entanglement, 
encoded by T,^^, is important. The complete cross section, ad + ae, is displayed in Fig. 3 (c). The entanglement 
contribution seems to change the cross section mainly along the rim of Fig. 3 (b). One sees the entanglement 
contribution clearer in Fig. 6 (e) where ad + ad is plotted together with ad- The closeness of the shape of |71p and 
2'!Tp^/^^ , plotted in Figs. 6 (f) and (b), is remarkable. 

Note that the entanglement contribution to the cross section, the second line of Eqs. (43), is positive deflnite in 
agreement of the notation, |7^P, used. Such a definite feature of the entanglement contribution appears natural since 
the asymptotic particle-holes states of the gas induced by the colliding particles should always increase the cross 
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section. 

Finally, the ratio {aj. + (Te)/<Jd, depicted in Fig. 7, shows the relative importance of the entanglement contributions 
around the rim of (7^, where the particle-hole components of the exchanged photon are around the Fermi surface. 
The nontrivial final states of the environment, taking into account by the entanglement contribution, represent the 
dominant contribution to the cross section in this kinematical region. 



D. Classicality 

There are two known necessary conditions of quantum-classical crossover, irreversibility and decoherence. As 
mentioned above, the latter corresponds to a property of the state at a given time and its building up in time can be 
seen as consistency. 

To simplify matter, let us consider n nonrclativistic charges. The generator functional like the one defined by Eq. 
(17) for the reduced density matrix of the charges is 

_D[A]D[a;]e5'^'"^'^^''^'^^'^'^^l"*'^["' '^'^•^''^^'^'"'l'*'^^^^ (84) 

where the boundary condition in time is OTP for particle trajectories x and CTP for the environment, the elec- 
tromagnetic field, and the electron gas. The electric current density of the given particle trajectories is denoted by 
j'^iu) = j'^bj'T^] ^-nd the external sources are suppressed for simplicity. Integration over the environment variables 
yields the efi^ective theory [36] 

Z = J D[x]e'^''ff^'^^ (85) 

where the effective action, 

Seff[x] = S[x+] - S[x-]+W''[-ej[x+],ej[x-]], (86) 
contains the influence functional, the generator functional of the electromagnetic field, 

Wij] = / dxdyf{x)D^,{x, y)f{y) + 0(?). (87) 
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The form (29) for the photon propagator yields for the quadratic part 

sftl^^(2)[3] ^ _^ Jdxdy[j'^{x)D;^{x,y)r{y)-r{x)D;,ix,y)f{y)] 
cj^7(2)[3-] = _2e^l dxdyr{x)Dl,{x,y)r{y) 



(88) 



where the parametrization = j/2 ± j is used. It is remarkable that the imaginary part is positive semidefinite, 
cjH^7(2)[j] > 0. The actual form is 



<^W'^[-ej[x+],ej[x-]] = -2e 



I 



dq 
(2^ 



i.nj{q)f+Dl{q){\j{q)f-\j{q)nf) 



(89) 



where the current in the right-hand side is j^{y) = j[y; a;+] — j[y; x~]. 

The functions Df and arc shown in Figs. 8. The sharp rim of D}. at the same position as in Fig. 7 indicates that 
the quasiparticles contributing strongly to entanglement scattering make the colliding particle trajectories consistent. 
The longitudinal, Coulomb-like contribution to consistency, shown in Fig. 8 (a) , is weak at small energy-momentum 
transfer due to screening. In this regime the transverse radiation field starts to generate consistency according to 
Fig. 8 (b). One expects an enhancement of these effects by further radiative corrections not considered in this work, 
including soft photons in the final state. 

Irreversibility is generated by the finite life time of quasi particles which is inversely proportional to the imaginary 
part of the inverse of the Feynman propagator at the quasiparticle dispersion relation. The CTP structure (29) of the 
photon propagator assures that both consistency and irreversibility are governed by the same dynamics, comprised 
in D*. The main lesson is that scattering processes with large values of D^, having fast-decaying quasiparticles and 
thereby strong irreversibility, tend to be classical. 



VII. SUMMARY 



A collision process with open environment channels in the asymptotic out-state sector, in particular, electron- 
proton collision in a nondegenerate electron gas, is treated in this work within the framework of the CTP formalism. 
The transition probability is obtained by the simple repetition of steps, followed in the reduction formulas and the 
environment is taken into account by using Schwinger-Dyson resummed photon propagator containing the one-loop 
self-energy. The dynamics of asymptotic environment states arises from the back-reaction of the collision process on 
its environment and can easily be handled by the algebraic solution of the Schwinger-Dyson-Kadanoff-Baym equation. 

It is found that back-reaction is important and dominates the cross section when the exchanged energy is around 
the Fermi level. Asymptotic environment states and back-reaction yield strong colliding particle-environment entan- 
glement in this regime. Hence, whenever soft back-reaction is dominant the collision is closer to being classical and is 
irreversible. 

One may wonder how further partial rcsummation of the perturbation series changes the results. Higher-order 
terms in the photon self-energy include multiple particle-hole pairs and, therefore, should enhance entanglement and 
classical features. It remains to be seen how vertex corrections modify the results. It seems reasonable to expect 
that our results remain qualitatively similar for degenerate electron gas and finite temperature effects should further 
strengthen entanglement and classical behavior. 

We believe that the results make the reevaluation of cross sections necessary for elementary processes in laboratory 
and in extraterrestrial plasma if their kinematical regime is close to the scale of the environment. Furthermore, this 
method may ultimately lead to an improved description of collision processes where the multiparticle aspects of the 
beam and the target can be taken into account and the construction of more powerful phenomenological kinetic models 
to describe the quantum-classical crossover. 
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